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INTRODUCTION 


Random  Number  Generators  (RNG 


~Tive  external  subroutines  are  available  to  generate  uncorrelated, 
floating  point  random  numbers  which  are  uniformly  distributed  between 
0.0  and  1.0.  The  subroutines  are  extremely  fast  (8  per  word)  and 
do  not  need  an  external  starting  number.  Additional  4ntry  points  are 
provided  so  that  a generated  sequence  of  randan  numbers  may  be  continued 
from  program  to  program.  | 

Procedure  for  RNG  - — ' 


The  Random  Number  Generator  subroutines  are  called  as  follows 


This  will  result  in  the  return  of  a floating  point  decimal  random 
number  between  0.0  and  1.0  to  core  location  R.  (See  example).  Each 
new  call  will  provide  a new  number. 


*A  large  part  of  the  material  for  the  assembly  language  generators 
cones  from  a Yale  Computer  Center  Memo.  No.  27G  by  J.  Lach,  and  fran 
an  IEEE  Spectrum  Article  in  Feb.  1967  by  R.  P.  Chambers,  "Randan  Number 
Generation  in  Digital  Computers”. 
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The  statements: 

CALL  GAS2  (K) 

CALL  GAS5  (R) 

may  be  used  to  extract  numbers  from  the  other  generators  in  the  same 
fashion.  Each  generator  has  a different  starting  number  and  constant 
with  the  result  that  all  five  generators  are  uncorrelated. 

In  order  to  continue  the  sequence  in  an  unbroken  fashion  in  a 
subsequent  program,  two  additional  entry  points  are  supplied  with 
each  generator. 

The  statements : 

CALL  SGAS1  (K) 

will  place  the  last  random  number  generated  as  a fixed  point  number 
in  core  location  R.  This  number  may  be  printed  or  punched  to  restart 
the  generator.  To  insure  that  all  significant  portions  are  transferred, 
a FORMAT  designation  of  012  should  be  used.  If  the  number  is  read  into 
the  next  program  In  the  same  format,  then  the  statement: 

CALL  U1AS1  (R) 

will  store  the  number  in  the  GAS1  generator  so  subsequent  calls  on 
GAS1  will  keep  the  sequence  unbroken  (see  example). 

Spec lal  Notes : 

1.  Therandom  number  generators  are  of  the  multiplicative,  congruential 
type.  * This  technique  employs  a recursion  relation  of  the  form: 

Kn+1  * r (Hn) 

where  Rn+i  is  the  (nd)1*1  rnndcro  number  generated  end  Y is  a function 
of  the  previous  random  number,  Hn.  In  particular,  the  function 
programmed  is 

RIU1  =■  C Kn  (Modulo 
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where  C ~ \j  . ^ and  obeys  the  condition 

C ■=  8 m ♦_  i 

1 2 

where  m is  an  Integer. 


It  should  be  kept  in  mind  that  true  randan  numbers  can  only  be 
generated  by  « random  process.  The  numbers  generated  by  a computer 
are  only  approximations  to  such  a process.  The  hope  is  that  the 
approximations  are  close  enough  for  the  calculations  at  hand. 


j.  Two  important  tesij  for  a random  number  generator  are  the  tests 
for  uniformity  of  the  distribution  between  0 and  1 and  for  the 
correlation  between  samples.  The  generators  used  here  have  been 
tested  for  correlation’  with  1000  samples  and  have  been  within  the 
expected  statistical  accuracy.  The  test  for  uniformity  has  been  made 
indir-’tly  by  a chi-square  test  for  normality1*  on  samples  produced 
by  this  generator  but  transformed  by  a "Statpack"*5  routine,  to  a 
normal  distribution.  For  500  samples  the  probability  of  exceeding 
chi-square  was  2%  at  a 9511  confidence  level.  The  generator  included 
in  the  Math-Pack  ,fc>  library  could  not  do  better  than  under  the 
same  conditions. 


4.  In  order  to  obtain  normal  distributions  this  generator  should  be 
used  in  conjunction  with  an  inverse  normal  distribution  such  as 
"•tatpack"^  provides.  (See  example). 

5.  Bach  generator  uses  a total  of  19  core  locations. 

6.  The  generators  are  in  File  3 of  CUR  Tape  UI83. 


A.  The  generation  of  100  uniformly  distributed  numbers  with  the 
last  number  saved  for  Program  a. 


DIMENSION  X(100) 
DO  101  = 1,  100 

10  CALL  (IASI  (X(l)) 
CALL  SGA31  (FIN) 
PUNCH  1,  FIN 
1 FORMAT  (012) 


(s’  Declare  Array  X 

Fill  X Array 
«'  Save  last  No. 

W Punch  last  No. 
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B.  The  generation  of  lOOnormally  distributed  randan  numbers  with 
zero  mean  and  unit  variance  starting  from  the  last  no.  of  Program  A. 

DIMENSION  X (100) 

1 FORMAT  (012) 

READ  (3,1)  START  0 Read  starting  No. 

CALL  LGAS1  (START)  & Enter  starting  No. 

DO  101  = 1,  100 

( 2 CALL  GAS1  (x(l))  Obtain  uniform  Nos. 

10  X(I)  = TINORM6  (X(I),  $2)  <sl  Obtain  normally  dist.  No. 

u 

N 

Conclusions 


The  randan  number  generators  presented  are  extremely  fast  and 
have  excellent  statistical  properties.  By  using  several  of  the 
generators  in  the  same  program  even  better  properties  may  be  obtained. 
The  main  use  will  probably  be  in  Monte  Carlo  Simulations. 


I 

i 


References 

1.  Chambers,  R.  P. , "Random-Number  Generation  on  Digital  Computers", 
IEEE  Spectrum,  Feb.  1967,  pp.  48-56. 

2.  Lach,  J. , "YCC  Memorandum  2 G",  Oct  1,  1964. 

3.  Newman,  H. , private  correspondence. 

4.  Foreman,  L. , private  correspondence. 

5.  Bauer,  W.  F.,  "The  Monte  Carlo  Method”,  J.  Soc.  Ind.  Appl.  Math., 
Vol.  6,  pp.  438-451,  1958. 

6.  Uni vac  1108  Math-Pack  and  Stat-Pack  Subroutines  UP-7542,  UP-7502. 


4 


‘ 


Cl  ilSIMCtriom 


NUMBER  or  CORIES  MADE 


UNCLASSIFIED 


|A,  Nuttall . Code  2020 


J.  Alflcro,  Code  20' 


J.  Almeida.  Code  2 


R.  Arnold.  Code  2070 


w.  As)iinWwici*  Coie  .2072 


207? 


JU.A*  Breslln,  Code  2072 


Cable.  Cole  2071 


R.  D.  Clark.  Code  207? 


W.  R.  Cleary  (UNIVAC) 


IE.  M.  Podeszwa.  Code  ?040 


D.  W«  Counsellor.  Code  2072 


R.  D.  Drinkarl.  Code  207? 


L.  T.  r In at e i n.  Code 


